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Abstract 



A covariant scalar-tensor-vector gravity theory is developed which allows 
the gravitational constant G, a vector field coupling u and the vector field mass 
\i to vary with space and time. The equations of motion for a test particle 
lead to a modified gravitational acceleration law that can fit galaxy rotation 
curves and cluster data without non-baryonic dark matter. The theory is 
' consistent with solar system observational tests. The linear evolutions of the 

metric, vector field and scalar field perturbations and their consequences for 
the observations of the cosmic microwave background are investigated. 
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99; 1 Introduction 

> ■ 

^ ■ Two theories of gravity called the nonsymmetric gravity theory (NGT) and the 

metric-skew-tensor gravity (MSTG) theory [2] have been proposed to explain the 
rotational velocity curves of galaxies, clusters of galaxies and cosmology without 
dark matter. A fitting routine for galaxy rotation curves has been used to fit a large 
number of galaxy rotational velocity curve data, including low surface brightness 
(LSB), high surface brightness (HSB) and dwarf galaxies with both photometric data 
and a two-parameter core model without non-baryonic dark matter [21 EJ. The fits to 
the data are remarkably good and for the photometric data only the one parameter, 
the mass-to-light ratio (M/ L) , is used for the fitting, once two parameters M and 
r are universally fixed for galaxies and dwarf galaxies. The fits are close to those 
obtained from Milgrom's MOND acceleration law in all cases considered [3]. A large 
sample of X-ray mass profile cluster data has also been fitted |HJ . 

The gravity theories require that Newton's constant G, the coupling constant 
7 C that measures the strength of the coupling of the skew field to matter and the 
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mass fj, of the skew field, vary with distance and time, so that agreement with the 
solar system and the binary pulsar PSR 1913+16 data can be achieved, as well as 
fits to galaxy rotation curve data and galaxy cluster data. In ref. the variation 
of these constants was based on a renormalization group (RG) flow description of 
quantum gravity theory formulated in terms of an effective classical action jH] . Large 
infrared renormalization effects can cause the effective G, j c , \i and the cosmological 
constant A to run with momentum k and a cutoff procedure leads to a space and 
time varying G, 7 C and /i, where \x = 1/ro and tq is the effective range of the skew 
symmetric field. 

In the following, we shall pursue an alternative relativistic gravity theory based 
on scalar-tensor- vector gravity (STVG), in which G, a vector field coupling constant 
uj and the mass fi of the vector field are dynamical scalar fields that allow for an 
effective description of the variation of these "constants" with space and time. We 
shall not presently consider the variation of the cosmological constant A with space 
and time. 

The gravity theory leads to the same modified acceleration law obtained from 
NGT and MSTG for weak gravitational fields and the same fits to galaxy rotation 
curve and galaxy cluster data, as well as to agreement with the solar system and 
pulsar PSR 1913+16 observations. An important constraint on gravity theories is 
the bounds obtained from weak equivalence principle tests and the existence of a 
"fifth" force, due to the exchange of a massive vector boson [7]. These bounds 
are only useful for distances < 100 A. U. and they cannot rule out gravity theories 
that violate the weak equivalence principle or contain a fifth force at galactic and 
cosmological distance scales. Since the variation of G in our modified gravity theory 
leads to consistency with solar system data, then we can explore the consequences 
of our STVG theory without violating any known local observational constraints. 

An important feature of the NGT, MSTG and STVG theories is that the modified 
acceleration law for weak gravitational fields has a repulsive Yukawa force added to 
the Newtonian acceleration law. This corresponds to the exchange of a massive 
spin 1 boson, whose effective mass and coupling to matter can vary with distance 
scale. A scalar component added to the Newtonian force law would correspond to an 
attractive Yukawa force and the exchange of a spin particle. The latter acceleration 
law cannot lead to a satisfactory fit to galaxy rotation curves and galaxy cluster data. 

In Section 8, we investigate a cosmological solution based on a homogeneous and 
isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime. We present 
a solution that can possibly fit the acoustic peaks in the CMB power spectrum 
by avoiding significant suppression of the baryon perturbations jH], and which can 
possibly be made to fit the recent combined satellite data for the power spectrum 
without non-baryonic dark matter. 

All the current applications of the three gravity theories that can be directly 
confronted with experiment are based on weak gravitational fields. To distinguish 
the theories, it will be necessary to obtain experimental data for strong gravita- 
tional fields e.g. black holes. Moreover, confronting the theories with cosmological 
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data may also allow a falsification of the gravity theories. Recently, the NGT and 
MSTG were studied to derive quantum fluctuations in the early universe from an 
inflationary-type scenario jUJ- 



2 Action and Field Equations 

Our action takes the form 

S = Scrav + S<h + Ss + Sm, 



where 
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Here, we have chosen units with c = 1, V M denotes the covariant derivative with 
respect to the metric g^ v . We adopt the metric signature r\^ v = diag(l, —1,-1,-1) 
where rj^ v is the Minkowski spacetime metric. We have 



Pj pA _|_ pA per per pA 

lA/J- /iA T 1 | t 2,l Act 1 Ul j, 



where denotes the Christoffel connection: 



^g Xa (d u g^ + d^gva - d a g^), 



(5) 



(6) 



and R = g^ u R^ u . Moreover, V(<f>) denotes a potential for the vector field while 
V(G), V{oj) and V(/x) denote the three potentials associated with the three scalar 
fields G(x),lu(x) and respectively. The field u(x) is dimensionless and A 

denotes the cosmological constant. Moreover, 
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The total energy-momentum tensor is given by 



T, 



/IV 



T. 



Mfiu 



S/ivi 



where T MfMU and T^ v denote the ordinary matter energy-momentum tensor and 
the energy-momentum tensor contribution of the M field, respectively, while Ts^ u 
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denotes the scalar G, to and \i contributions to the energy-momentum tensor. We 
have 



2 5S 



M 



-Ti 



2 5Sa 



The matter current density J M is defined in terms of the matter action Sm '■ 

1 SS M _ _ JPu 



-T, 



2 5S, 



-Ti 



V~9 dPti 

We obtain from the variation of g^ u , the field equations 
where = R^ v — \g pv R. We have 

Q/iu = G(v a v a e 9llu - v M v,e), 



(9) 



(10) 



(12) 



where 9(x) = 1/G(x). The quantity Q^ v results from a boundary contribution 
arising from the presence of second derivatives of the metric tensor in R in Scrav 
These boundary contributions are equivalent to those that occur in Brans-Dicke 
gravity theory [TJIj. We also have 



T^^y — LO 



B«B va - 9flv [ -B p °B pa + V((j>)\ + 2 dV{<f)) 



dg 



The G[x) field yields the energy-momentum tensor: 
Tg^u '- 
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V M GV„G - 2^^ - g» u (^V a GV a G - V(G)j 



dg 



fJbV 



Similar expressions can be obtained for T wfiu and T^^. 
From the Bianchi identities 

V V G» V = 0, 



and from the field equations (fTTf . we obtain 



G o7rG 



V V Q^ = 0. 



A variation with respect to 0» yields the equations 



V V B» V + 



dV(</>) 



+ -V v toB 



LO 



LO 



(13) 
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(15) 



(16) 
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Taking the divergence of both sides with respect to V„, we get 



J 
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-V u ujB 



-J* 
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If we assume that the current density J M is conserved, we obtain the equation 

In standard Maxwell-Proca theory, the conservation of the current J M is a separate 
physical assumption. 

We shall choose for the potential V(0): 

V(0) = -^ 2 0^ + I^(0), (20) 

where W(4>) denotes a vector field M self- interaction contribution. We can choose 
as a model for the self-interaction: 

W(<j>) = \g{^)\ (21) 

where g is a coupling constant. 

The effective gravitational "constant" G(x) satisfies the field equations 



QV a GV a fi + -V a fiV a fi + uui/G^^- 

and the last term arises from the \i dependence of V(4>) in (J2UJ). 
If we adopt the condition 



V Q V a G + V'(G) + N= 1 -G 2 (t + J^J , 

where 

N = | (~ V a GV a G -V{G)^+G (U/ a uV a u - 

-3eV Q GV Q G + — ^ V a /iV> - V(jjl)) + Y^^" 9 ' ( 23 ) 
and T = g^ u T^ u . The scalar field cu(x) obeys the field equations 

V v Vu + V'(u) + F = 0, (24) 

where 

F = -0 V a GV a cu + G (-B^B^ + V(<P)^j . (25) 
The field /i(x) satisfies the equations 

V a V> + V'{n) + P = 0, (26) 

where 



(27) 



= 0, (28) 



then (fTTj) takes the form 



VV^ - R^fo + fi 2 <P, - - -VwB F = -J„. (29) 

The test particle action is given by 

/r dx^ 
dr — A J druuf)^—-—, (30) 

where r is the proper time along the world line of the test particle and m and A 
denote the test particle mass and coupling constant, respectively. The stationarity 
condition SStp/Sx^ = yields the equations of motion for the test particle 

where 

F = *t»B* v ^ + AV" W (^) - AV a u,(>^). (32) 

The action for the field B^ v is of the Maxwell-Proca form for a massive vector 
field 0^. It can be proved that this theory possesses a stable vacuum and the 
Hamiltonian is bounded from below. Even though the action is not gauge invariant, 
it can be shown that the longitudinal mode (fio (where M = (0o, <f>%) {i = 1,2,3)) 
does not propagate and the theory is free of ghosts. Similar arguments apply to the 
MSTG theory 0. 

The Hamilton-Dirac (HD) method is a tool for investigating the constraints and 
the degrees of freedom of a field theory [TI]. The HD procedure checks the theory 
for consistency by producing the explicit constraints, and counting the number of 
degrees of freedom. It is a canonical initial value analysis. When the field theory 
is coupled dynamically to gravity, many vector field theories are ruled out, because 
the constraints produced by the canonical, Cauchy initial value formalism yield 
"derivative coupled" theories (the Christoffel connections do not cancel out). The 
Maxwell and Maxwell-Proca theories are prime examples ^2] of consistent vector 
field theories. With or without the gravitational field coupling, Maxwell's theory 
has two degrees of freedom and Maxwell-Proca has three, and they are stable and 
satisfy a consistent Cauchy evolution analysis. Many other vector field theories are 
derivative coupled when the gravitational field is introduced into the action. Severe 
singularity problems can appear in vector-gravity coupled theories, which render 
them inconsistent. There are no pathological singularities in the Maxwell-Proca 
theory coupled to gravity, when one solves for the second time derivative in the 
canonical initial value formulation. In other vector theories, singularities occur that 
spoil the stability of the theory and rule them out as physically unviable theories. 

It is possible to attribute the mass /i to a spontaneous symmetry breaking mech- 
anism, but we shall not pursue this possibility at present. 
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3 Equations of Motion, Weak Fields and the Mod- 
ified Gravitational Acceleration 

Let us assume that we are in a distance scale regime for which the fields G, uj and 
\i take their approximate renormalized constant values: 

G ~ Gq(1 + Z), uj ~ ujqA, /i ~ /io-B, (33) 

where G , uj q and /x denote the "bare" values of G, u and /x, respectively, and Z, A 
and B are the associated renormalization constants. 

For a static spherically symmetric field the line element is given by 

ds 2 = 7(r)rft 2 - a(r)dr 2 - r 2 (d6 2 + sin 2 6d(f) 2 ). (34) 

The equations of motion for a test particle obtained from (jMJ), and (|3l3|) are 
given by 

d 2 r a' fdr\ 2 rfd8\ 2 f sin 2 9 \ f ad(j) \ 2 i f dt\ 2 



dT 2 + 2a\dTj aXdr) V a )\ dr J + 2a\dr) 
where a = Xu/m. 

The orbit of the test particle can be shown to lie in a plane and by an appropriate 
choice of axes, we can make 6 = ir/2. Integrating Eq.(|38|) gives 

r — = J, (39) 
dr 

where J is the conserved orbital angular momentum. Integration of Eq. gives 

^ = --M> + £), (40) 
dr 7 

where E is the constant energy per unit mass. 

By substituting (j4*Uj) into and using (|3*9*jl. we obtain 

cPr a' fdr\ 2 J 2 j' , 1 /d0 o \, , 



We do not have an exact, spherically symmetric static solution to our field equa- 
tions for a non-zero V{4>). However, if we neglect V(4>) in Eq. (jl7j) and A in Eq.(|llj). 
then the exact static, spherically symmetric (Reissner-Nordstrom) solution in empty 
space yields the line element 



, . 2GM Q 2 \ , ., ( 2GM Q 2 \ 1 , 
ds 2 =[l — + -^)dt 2 -(l h-^-J dr 1 -r*(d0* + sin 2 6d^), (42) 



where M is a constant of integration and 

Q 2 = AitGuje 2 . (43) 
Here, e denotes the "charge" of the spin -1 vector particle given by 

d 3 xq, (44) 



where the matter current density J M is identified as J M = (q, J % ). 

For large enough values of r, the solution (J42j) approximates the Schwarzschild 
metric components a and 7: 

1 . , 2GM . . 

a(r) ~l-2GM/r' 7(0-1- — (45) 

It is not unreasonable to expect that a static, spherically symmetric solution of the 
field equations including the mass term /1 in the field equations (|17|) will approximate 
for large values of r the Schwarzschild metric components (|45|). 

We assume that 2GM/r <C 1 and the slow motion approximation dr/ds ~ 
dr/dt <C 1. Then for material test particles, we obtain from (J3Tj) . (jUJ) 
and (EE5J) : 

cPr J 2 GM d(j)o 

^-^ + — = ^> ^ 

where Jn is the Newtonian orbital angular momentum. 

For weak gravitational fields to first order, the static equations for 0o obtained 
from (|29j) are given for the source-free case by 

V 2 o - v 2 <po = 0, (47) 

where V 2 is the Laplacian operator, and we have neglected any contribution from 
the self-interaction potential W{4>). For a spherically symmetric static field <fto, we 
obtain 



% + V - tffo = 0. (4* 



This has the Yukawa solution 



Mr) = -£*tel t (49) 
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where (3 is a constant. We obtain from (|46j) : 



rfV _ 4 + GM = ^ exp(-//r) (1 + }) 
c?t 2 r 3 r 2 r 2 



(50) 



where = a (3. 

We observe that the additional Yukawa force term in Eq. (j50|) is repulsive in 
accordance with the exchange of a spin 1 massive boson. We shall find that this 
repulsive component of the gravitational field is necessary to obtain a fit to galaxy 
rotation curves. 

We shall write for the radial acceleration derived from f!5()[) : 



a{r) 



+ g exp(-r/r )/ 1+ r\ 
r 2 \ r J 



and Goo is defined to be the effective gravitational constant at infinity 



(51) 



(52) 



Here, M denotes a parameter that vanishes when u = and Go is Newton's gravi- 
tational "bare" constant. The constant K is chosen to be 



K = G \ MM . 



(53) 



The choice of K, which determines the strength of the coupling of B^ u to matter 
and the magnitude of the Yukawa force modification of weak Newtonian gravity, 
is based on phenomenology and is not at present derivable from the STVG action 
formalism. 

By using (J52j) . we can rewrite the acceleration in the form 



a(r) 



G M 



'Mo 
M 



1 - exp(-r/r ) H 

r 



(54) 



We can generalize this to the case of a mass distribution by replacing the factor 
GqM/t 2 in (|54j) by GoAi(r)/r 2 . The rotational velocity of a star v c is obtained from 
v 2 (r)/r = a(r) and is given by 



'GnM(r) 



i Mo 
M 



1 - exp(-r/r ) ( 1 + 



r 



1/2 



(55) 



The gravitational potential for a point source obtained from the modified accelera- 
tion law (|54|) is given by 



$(r) 



G M 



1 + 



I Mo 
M 



1 - exp(-r/r )) 



(56) 



The acceleration law (Jo4J) can be written as 

G(r)M 



a(r) 



where 



G(r) = G c 



1 + 



'Mo 
M 



exp(-r/r ) ^1 + 



(57) 



(5* 



is an effective expression for the variation of G with respect to r. A good fit to a 
large number of galaxies has been achieved with the parameters E] : 



9.60 x 10 11 M, 



r = 13.92 kpc = 4.30 x 10 22 cm 



0; 



(59) 



In the fitting of the galaxy rotation curves for both LSB and HSB galaxies, using 
photometric data to determine the mass distribution Ai(r) 3J, only the mass-to- 
light ratio (M/ L) is employed, once the values of M and r are fixed universally for 
all LSB and HSB galaxies. Dwarf galaxies are also fitted with the parameters j^j: 



M = 2.40 x 10 11 M Q , r = 6.96 kpc = 2.15 x 10™ cm. 



22 



(60) 



By choosing values for the parameters Goo, (Mo) c iust and (r ) c iust, we are able to 
obtain satisfactory fits to a large sample of X-ray cluster data j^j. 



4 Orbital Equations of Motion 

We set 9 = n/2 in divide the resulting expression by dr 2 and use Eqs.fl3^|) and 
(l4Tty to obtain 

fdr\ 2 J 2 1 . , E , . 

— + — (700 + E 2 = . 61 

var / ar z cey a 

We have ds 2 = Edr 2 , so that ds/dr is a constant. For material particles E > and 
for massless photons E = 0. 

Let us set u = 1/r and by using (JHHJ), we have dr/dr = —Jdu/dcp. Substituting 
this into (|6ip. we obtain 

On substituting ()45|) and dr/dcj) = —(l/u 2 )du/d<j) into (p)2*j) . we get after some ma- 
nipulation: 

d 2 u EGM EK ( 1 \/ 1 \ 2 
— +m=— 7r ex P 1 + ) + 3GMu , (63) 



J 2 J 2 \ TqUJ \ TqU 

where r = 



10 



For material test particles E — 1 and we obtain 

d 2 u GM nm 2 K ( 1 \/ 1 \ 

1- 7/ = 1- 3G A/f?/ 2 p^ri 1-1 



+ u= — - + 3GMu 2 - — exp 1 + . (64) 

ui/j- J z J z V tquJ \ r uj 

On the other hand, for massless photons ds 2 = and E = and (JfiHj) gives 

d 2 v 

— + it = 3GMu 2 . (65) 

5 Solar System and Binary Pulsar Observations 

We obtain from Eq. (}64|) the orbit equation (we reinsert the speed of light c): 



d 2 u GM K . . . 

^ + M = ^~^ eXp( - r/ro) 



l+(^ 



Using the large r weak field approximation, and the expansion 



3GM 2 
H 5— w • (66) 



exp(-r/r ) = 1 - — + \(— ) + ... (67) 



r 2\r 

we obtain the orbit equation for r <^ r : 

d 2 u „ GM 



(i0 2 ' c 2 
where 

GM K , . 

^ = TPf - "T72-- ( 69 ) 

We can solve Eq. (J68|) by perturbation theory and find for the perihelion advance 
of a planetary orbit 

Au=^-(GM Q -K Q ), (70) 

c ±j 

where Jjy = (GMqL/c 2 ) 1 ^ 2 , L = a(l — e 2 ) and a and e denote the semimajor axis 
and the eccentricity of the planetary orbit, respectively. 

We now use the running of the effective gravitational coupling constant G = 
G(r), determined by (p)%j) and find that for the solar system r <C r , we have G ~ Go 
within the experimental errors for the measurement of Newton's constant Go- We 
choose for the solar system 

1.5 km (71) 

cr 

and use G = Go to obtain from (J7U)) a perihelion advance of Mercury in agreement 
with GR. The bound (J71)) requires that the coupling constant u varies with distance 
in such a way that it is sufficiently small in the solar system regime and determines 
a value for M , in Eq. fjo^j) . that is in accord with the bound (J7TJ). 
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For terrestrial experiments and orbits of satellites, we have also that G ~ Go and 
for sufficiently small, we then achieve agreement with all gravitational terrestrial 
experiments including Eotvos free-fall experiments and "fifth force" experiments. 

For the binary pulsar PSR 1913+16 the formula (fTUj) can be adapted to the 
periastron shift of a binary system. Combining this with the STVG gravitational 
wave radiation formula, which will approximate closely the GR formula, we can 
obtain agreement with the observations for the binary pulsar. The mean orbital 
radius for the binary pulsar is equal to the projected semi-major axis of the binary, 
(r) N = 7 x 10 10 cm, and we choose {r) N <C r . Thus, for G = G within the exper- 
imental errors, we obtain agreement with the binary pulsar data for the periastron 
shift when 

4.2 km. (72) 



For a massless photon E = and we have 

d 2 u GM 



+ u = 3—^ u \ (73) 



dcf) 2 c 

For the solar system using G = Go within the experimental errors gives the light 
deflection: 

A S = ^ (74, 

in agreement with GR. 



6 Galaxy Clusters and Lensing 

The bending angle of a light ray as it passes near a massive system along an ap- 
proximately straight path is given to lowest order in v 2 /c 2 by 

= ^ J \a L \dz, (75) 

where _L denotes the perpendicular component to the ray's direction, and dz is the 
element of length along the ray and a denotes the acceleration. 
From (JTSJ), we obtain the light deflection 

4GM 4G M 

A = ^ = ^iT' (76) 

where 



M = M(1 + a/^). (77) 



M 

The value of M follows from (|HH|) for clusters as r > ro and 



G(r)^G oo = G (l + J^). 178) 
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We choose for a cluster M = 3.6 x 1O 15 M and a cluster mass M c i ust ~ 1O 14 M , 
and obtain 

M ) dust ^ ^ 

We see that M ~ 7M and we can explain the increase in the light bending without 
exotic dark matter. 

From the formula Eq. ()54)1 for r > ro we get 

( \ GoM (art 
a ( r ) = —■ ( 8 °) 

We expect to obtain from this result a satisfactory description of lensing phenomena 
using Eq. ffTo'j) . 

An analysis of a large number of clusters shows that the MSTG and STVG 
theories fit well the cluster data in terms of the cluster mass, M clust , and an average 
value for the parameter Mq j^J. 



7 Running of the Effective Constants G, uj and fi 

The scaling with distance of the effective gravitational constant G, the effective 
coupling constant uj and the effective mass \i is seen to play an important role in de- 
scribing consistently the solar system and the galaxy and cluster dynamics, without 
the postulate of exotic dark matter. We have to solve the field equations J22J), (jHI 
and (|2l)j) with given potentials V(G),V(oj) and V(/x) to determine the variation of 
the effective constants with space and time. These equations are complicated, so 
we shall make simplifying approximations. In Eq. (J22j) . we shall neglect the contri- 
butions from N and obtain for T = A = 0: 

V U VG + V\G) = 0, (81) 

where f'(y) = df/dy. The effective variation of G with r is determined by Eq.(|58j). 
We obtain from (JHTj) for the static spherically symmetric equations 

V 2 G(r) - V'(G) = G"(r) + -G'ir) - V'(G) = 0. (82) 

r 

By choosing the potential 

^(SXSW-^K-S- (83) 

we obtain a solution to (|82p for G(r) given by (|58|). The neglect of the contribu- 
tions N can only be justified by solving the complete set of coupled equations by a 
perturbation calculation. We will not attempt to do this in the present work, but 
we plan to investigate this issue in a future publication. 
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We see from fjHHjl that for r < r we obtain G{r) ~ Go- As the distance scale 
approaches the regime of the solar system r < 100A£7. where 1A.U. = 1.496 x 
10 13 cm = 4.85 x 10 _9 kpc, then (J54j) becomes the Newtonian acceleration law: 

in agreement with solar physics observations for the inner planets. 

Let us make the approximation of neglecting F((p) in Eq.(|24|). In the static 
spherically symmetric case this gives 

oj"(r) + -v'(r) - V'iuS) = 0. (85) 
r 

We choose as a solution for uj(r): 

uo(r) = uj {1 + ZJ[1 - exp(-/Zr)(l + ~pr)]}, (86) 
where ZJ and JZ are positive constants. The potential V(u) has the form 

V(u) = -^ 2 /Z 2 cJ 2 exp(-277r)(l + 2/Zr - JI 2 r 2 ). (87) 

For the variation of the renormalized mass \i = /i(r), we find that a satisfactory 
solution to Eq.(J21J) should correspond to a a /x(r) = l/r (r) that decreases from a 
value for the inner planets of the solar system, consistent with solar system obser- 
vations, to a small value corresponding to r for the galaxy fits, r = 14 kpc, and to 
an even smaller value for the cluster data fits. 

The spatial variations of G(r), u{r) and /i(r) = l/r (r) can be determined 
numerically from the equations (j22j), ([HI and (pjj) with given potentials V(G), V(u) 
and V(fj), such that for the solar system and the binary pulsar PSR 1913+16 the 
bounds (fTTf and (fT2|) are satisfied by the solutions for G(r),ui(r) and fi(r). The 
spatial variations of G, the coupling constant uj and the range vq are required to 
guarantee consistency with solar system observations. On the other hand, their 
increase at galactic and cosmological distance and time scales can account for galaxy 
rotation curves, cluster lensing and cosmology without non-baryonic dark matter. 

We have constructed a classical action for gravity that can be considered as an 
effective field theory description of an RG flow quantum gravity scenario as described 
in refs. jU] and [2|. 

The fitting of the solar system, galaxy and the clusters of galaxies data depends 
on the running of the of the "constants" G, u and tq. They should increase from 
one distance scale to the next according to the renormalization group flow diagrams, 
or the solutions of the classical field equations in the present article. In a future 
article, the author plans to provide a more complete determination of the running 
of the constants. However, the present article describes the basic scenario and the 
ideas underlying the theory. 
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8 Cosmology 



Let us now consider a cosmological solution to our STVG theory. We adopt a 
homogeneous and isotropic FLRW background geometry with the line element 



ds 2 = dt 2 - a 2 (t) 



dr 2 



kr 2 



+ r 2 d& 



(88) 



where dVL 2 = d6 2 + sin 2 9d(j) 2 and k — 0, —1, +1 for a spatially flat, open and closed 
universe, respectively In this background spacetime, we have 0o = i> 7^ 0, (pi = 
and B^ u = 0. 



We define the energy-momentum tensor for a perfect fluid by 



= (p + p)u^u v 



19) 



where 
we have 



dx^/ds is the 4- velocity of a fluid element and g^u^u 1 ' = 1. Moreover, 



p = Pm + P<f> + ps, P = Pm +P& + PS, 



(90) 



where pi and pi denote the components of density and pressure associated with the 
matter, the </> M field and the scalar fields G, u and p, respectively 
The Friedmann equations take the form 



H 2 {t) + 



a 2 (t) 



*irG(t)p(t) dG A 
3 a G 3 ' 



d(t) 



47rG(t) . , \ / w 1{G G 2 2dG\ A 



a(t) 

where H(t) = d(t)/a(t). 

Let us make the simplifying approximation for equations 



(91) 



(92) 



g + 3Hg + v\g) 



-G0Q 2 



p - 3p + 



A 



AnG g 



(93) 



where Q(t) = G(t)/G . A solution for Q in terms of a given potential V{Q) and for 
given values of p and p can be obtained from ()93|) . 

The solution for Q must satisfy a constraint at the time of big bang nucleosyn- 
thesis [131 ■ The number of relativistic degrees of freedom is very sensitive to the 
cosmic expansion rate at 1 MeV. This can be used to constrain the time dependence 
of G. Recent measurements of the 4 He mass fraction and the deuterium abundance 
at 1 MeV leads to the constraint G(t) ~ Go- We impose the conditions Q(t) — > 1 
as t — > tsBN and Q(t) — » 1 + uj as t — > ts^s where t^BN and ^5 denote the times 
of the big bang nucleosynthesis and the surface of last scattering, respectively. A 
possible solution for Q can take the form 



G{t) = l + cu 



1 - exp(-i/T) 



1 + l 



(94) 



15 



where u and T are constants and u is a measure of the magnitude of the scalar field 
if). We have for t » T that Q -> 1 + u and for t « T that Q 1. We get 

^^apH/T), ^ = ^ex P H/T)(l-i). (95) 

It follows that Q(t) ^ for t » T, which allows us for a suitable choice of T to 
satisfy the experimental bound from the Cassini spacecraft measurements [Tij : 

\G/G\ ~10- 13 yr _1 - (96) 

A linear perturbation on the FLRW background will link the theory with obser- 
vations of anisotropics in the CMB as well as galaxy clustering on large scales. The 
basic fields are perturbed around the background spacetime (denoted for a quantity 
7 by 7). In the conformal metric with the time transformation drj = dt/a(t): 



ds 2 = a 2 (ri)(dr] 2 - dx 2 ), (97) 

the metric perturbations are in the conformal Newtonian gauge 

g 00 (x,t) = a 2 (t)(l + 2$(x,t)), 9ij (x,t) = a 2 (t)(l - 2$(x,t))5 ij , (98) 

where $ is the gravitational potential. The vector field perturbations are defined by 

<j)„{x,t) = a(t)(<^(t) + <%(£,£)), (99) 

where 4>%{t) = 0. Denoting by Xi t ne scalar fields Xi — G, X2 = ^ and X3 — Hi the 
scalar field perturbations are 

Xi{^t)=Xi{t)+Sxi{x,t). (100) 



The problem that gravity theories such as MSTG and STVG face in describ- 
ing cosmology with no cold dark matter (CDM) (non-baryonic dark matter) is the 
damping of perturbations during the recombination era. In a pure baryonic universe 
evolving according to Einstein's gravitational field equations, the coupling of baryons 
to photons during the recombination era will suffer Silk damping, causing the colli- 
sional propagation of radiation from overdense to underdense regions [HI EI] • I n the 
CDM model, the perturbations Scdm are undamped during recombination, because 
the CDM particles interact with gravity and only weakly with matter (photons). 
The Newtonian potential in the CDM model is approximately given by 

k 2 <5> ~ A7cG (p b 5 b + PcdmScdm), (101) 

where k 2 denotes the square of the wave number, p b and Pcdm denote the densities 
of baryons and cold dark matter, respectively, and <5; denotes the perturbation den- 
sity contrast for each component i of matter. If pcdm is sufficiently large, then Scdm 
will not be erased, whereas S b decreases during recombination [16 . In STVG, the 
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imperfect fluid plasma before recombination has two components: the dominant neu- 
tral vector field component that does not couple to photons and the baryon-photon 
component. The vector field component has zero pressure and zero shear viscos- 
ity, so the vector field perturbations are not Silk damped like the baryon-radiation 
perturbations, for the latter have non-vanishing pressure and shear viscosity. 
In STVG, we have 

a = ^% n, = ^, n s = 8 -^. (102) 

b 3H 2 ' * 3H 2 6 3H 2 K J 

We also have a possible contribution from massive neutrinos 

a, = d03) 

where p v denotes the density of neutrinos. We assume that the vector field density 
dominates. 

The Newtonian potential in our modified gravitational model becomes 

/c 2 $ ~ 47rG Icn [p b 5 b + p v 5 v + + p s 6 s ] } (104) 

where G ren is the renormalized value of the gravitational constant. Assuming that 
the density p^ is significant before and during recombination, we can consider fit- 
ting the acoustic peaks in the power spectrum in a spatially flat universe with the 
parameters 

tt = tt b + Q, v + fi^ + Q, v + A = 1. (105) 

The fitting of the acoustic peaks in the CMB power spectrum does not per- 
mit a too large value of Q\. Moreover, the neutrino contribution is constrained 
by the three neutrinos having a mass < 2 eV. We could choose, for example, 
Q b = 0.04, ft^p = 0.25, Q u = 0.01,Oa — 0.7 as a possible choice of parameters to 
fit the data. There are now new data from the balloon borne Boomerang CMB 
observations that together with other ground based observations and WMAP 
data determine more accurately the height of the third acoustic peak in the angular 
CMB power spectrum JH]- The ratio of the height of the first peak to the second 
peak determines the baryon content Q b ~ 0.04. The height of the third peak is 
determined by the amount of cold dark matter, and in the modified gravity theory 
by the possible amounts of scalar ip, G, uo and p contributions. In particular, the 
dominant neutral ip vector component perturbations will not be washed out before 
recombination. 

Can the effects of gravitational constant renormalization together with the pos- 
sible effects of the densities p^ and ps describe a universe which can reproduce 
the current galaxy and CMB observations? The answer to this problem will be 
addressed in a future publication. 
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9 Conclusions 



A modified gravity theory based an a D = 4 pseudo-Riemannian metric, a spin 1 
vector field and a corresponding second-rank skew field and dynamical scalar 
fields G, uj and /i, yields a static spherically symmetric gravitational field with 
an added Yukawa potential and with an effective coupling strength and distance 
range. This modified acceleration law leads to remarkably good fits to a large 
number of galaxies [3] and galaxy clusters [3] . The previously published gravitational 
theories NGT [T] and MSTG [2] yielded the same modified weak gravitational field 
acceleration law and, therefore, the same successful fits to galaxy and cluster data. 
The MSTG and STVG gravity theories can both be identified generically as metric- 
skew-tensor gravity theories, for they both describe gravity as a metric theory with 
an additional degree of freedom associated with a skew field coupling to matter. 
For MSTG, this is a third-rank skew field F^ u x, while for STVG the skew field is a 
second-rank tensor B^ u . However, MSTG is distinguished from STVG as being the 
weak field approximation to the nonsymmetric gravitational theory (NGT). 

An action Ss for the scalar fields G(x), uj{x) and fi(x) = 1/tq(x) and the field 
equations resulting from a variation of the action, 5Ss = 0, can be incorporated into 
the NGT and MSTG theory. The dynamical solutions for the scalar fields give an 
effective description of the running of the constants in an RG flow quantum gravity 
scenario, in which strong infrared renormalization effects and increasing large scale 
spatial values of G and u lead, together with the modified acceleration law, to 
a satisfactory description of galaxy rotation curves and cluster dynamics without 
non-baryonic dark matter. 

We have demonstrated that a cosmological model with the renormalized gravi- 
tational constant G ren and contributions from the scalar fields G(t), u){t) and 
can possibly lead to a satisfactory description of the distribution of galaxies and the 
CMB power spectrum in a baryon dominated universe. 

The neutral vector particle <fr does not couple to radiation and it has zero pres- 
sure p and zero shear viscosity. Since it dominates the period of recombination, 
its perturbations associated with the plasma fluid will not be washed out by Silk 
damping. In contrast to standard dark matter models, we should not search for 
new stable particles such as weakly interacting massive particles (WIMPS) or neu- 
tralinos, because the fifth force charge in STVG that is the source of the neutral 
vector field (skew field) is carried by the known stable baryons (and electrons and 
neutrinos). This new charge is the source of a fifth force skew field that modifies 
the gravitational field in the universe. 
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